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For electron transport in parallel-plane semiconducting structures, a model is developed that 
unifies ballistic and diffusive transport and thus generalizes the Drude model. The unified model 
is valid for arbitrary magnitude of the mean free path and arbitrary shape of the conduction band 
edge profile. Universal formulas are obtained for the current-voltage characteristic in the nondegen- 
erate case and for the zero-bias conductance in the degenerate case, which describe in a transparent 
manner the interplay of ballistic and diffusive transport. The semiclassical approach is adopted, but 
quantum corrections allowing for tunneling are included. Examples are considered, in particular 
the case of chains of grains in poly crystalline or microcrystalline semiconductors with grain size 
comparable to, or smaller than, the mean free path. Substantial deviations of the results of the 
unified model from those of the ballistic thermionic-emission model and of the drift-diffusion model 
are found. The formulation of the model is one-dimensional, but it is argued that its results should 
not differ substantially from those of a fully three-dimensional treatment. 

PACS number(s): 05.60.-k, 05.60. Cd, 72.10.-d, 72.20.-i 

I. INTRODUCTION 

Electron transport in semiconducting structures is ballistic if the mean free path is much larger than the character- 
istic dimensions of the sample, and it is diffusive if the mean free path is much smaller than these. In the first case, 
there is no impurity or lattice scattering, and the current is determined by the ballistic motion in the electric field 
in the second case, scattering predominates and is described within the drift-diffusion scheme HQ. In many 
cases of physical relevance, however, the mean free path is neither large nor small compared with the characteristic 
dimensions of the sample. Thus, formulas for the current-voltage characteristic have appeared in the literature which 
combine features of the two limiting types of transport mechanism for particular conduction band edge profiles, e.g., 
for a single barrier 

In the present paper, we consider electron transport in parallel-plane semiconducting structures, i.e., structures 
whose parameters vary in one direction only. We develop a one-dimensional transport model which is valid for any 
magnitude of the mean free path and any form of the band edge profile, and thus unifies the ballistic and diffusive 
transport mechanisms. It is based on the idea that the electrons move ballistically in the electric field over intervals 
with average length equal to a universal mean free path, after which they are thermalized into a state of local 
equilibrium characterized by a quasi- Fermi level (electrochemical potential). The length of the sample is made up of 
random configurations of such ballistic intervals. Averaging over these configurations results in a unified description 
of electron transport, in which purely ballistic and purely diffusive transport appear as limiting cases. We work within 
the semiclassical approach [ Oft, which allows a concise and transparent formulation. However, quantum tunneling 
("thermionic field emission" |13| , |l4| ]) is taken into account in WKB approximation. 

The description of transport in terms of ballistic motion over intervals of the average length of the mean free path 
with thermalization at the end is, of course, also the basis of the Drude model |l2|JH]| , |l(| and of the relaxation-time 
approximation of the Boltzmann equation. There, however, the further development makes use of the assumption that 
the mean free path is small compared to the characteristic dimensions of the sample, leading to a diffusive description 
of the transport. In contrast to this, such an assumption is not made in the present work, and the magnitude of 
the mean free path relative to the sample dimensions determines the relative importance of the ballistic and diffusive 
transport mechanisms. The unified description thus is a generalization of the Drude model; it is particularly relevant 
to polycrystalline and microcrystalline materials (in the following indiscriminately referred to as "polycrystallinc 
materials") when the grain size is comparable to, or smaller than, the mean free path. In this case, the grains must 
not be considered separately and the sample must be treated as a whole. 

Our approach is the semi-phenomenological one commonly used in the description of transport in semiconductor 
devices (cf., e.g., Ref. Q). Its prominent and useful feature is that its results are obtained in closed form, which 
allows one to analyze the physics of the transport process in a transparent way. We have been able to formulate 
our model in one dimension. For transport across parallel-plane structures as considered in the present work, a 
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fully three-dimensional treatment of the thermalization process is expected not to change the resulting formulas in a 
substantial way (cf. Sec. III.D below for a more detailed discussion of this point). 

In the following section, we introduce the basic assumptions of the unified model and explain the procedure of 
averaging over the random configurations of ballistic intervals across which the electrons travel without scattering. 
In Sec. Ill, we present our principal result: a universal formula for the current-voltage characteristic reflecting the 
interplay of ballistic and diffusive transport. This formula generalizes previous expressions proposed for a single 
barrier. Its distinctive new feature is the appearance of a shape term that depends on the detailed structure of the 
band edge profile and also explicitly on the mean free path. This term is essential, e.g., for the description of transport 
in polycrystalline materials. Section IV deals with two numerical examples. First, we investigate chains of identical 
grains in polycrystalline materials. Here, the effect of the relative magnitudes of mean free path and characteristic 
length of the sample, i.e., the relative importance of the ballistic and diffusive mechanisms, is demonstrated explicitly. 
It is found that the (zero-bias) conductivity, defined as conductance times sample length, generally depends upon 
the number of grains in the sample. Second, as an example of a degenerate system, we discuss the grain barrier 
conductance of a highly doped polycrystalline material as a function of temperature. Section V contains a summary 
and some concluding remarks. 



II. THE MODEL 
A. Basic formulation 

The unified model is based on the following one-dimensional scheme (cf. Fig. |l|). The length of the sample (extending 
from x = to x — S) is covered by a chain of N intervals {i} with end points xf_ x and xf , i — 1,...,N (xf > 
x i-i't x o — 0: x n = across which the electrons move ballistically in the field of the conduction band edge potential 
E c (x). There may be any number of intervals in a chain, N — 1, oo. At the end points of these "ballistic intervals" 
(in particular, at the fixed end points of the sample), the electrons are equilibrated (thermalized) . These points of 
local equilibrium are theoretical constructs which may be thought to be connected via thin ideal leads (flat potentials) 
to large ideal reservoirs in equilibrium, with chemical potentials equal to the quasi- Fermi level at the points plf] . The 
current flow in a ballistic interval is assumed to result from the injection of electrons at the end points and their 
transmission through the interval, in line with Landauer's view of conduction as a transmission phenomenon [^8[-pl|. 




FIG. 1. Averaging over the ballistic configurations: single peak of the band edge profile E c [x) at x — X. 



The current in the ballistic interval {i} is given by 

/>oo 

j(i)=-e de [Tt(e + E c (x^ 1 ))f(x?_ 1 ;e)-T?i6 + E c (x l ))f(x?;e)] , i = l,...,N, (1) 

Jo 

where f(x; e) is the phase space density of the electrons at position x with kinetic energy e of the motion in the 
x-direction, and T^(E) is the (classical or quantal) probability for ballistic transmission at total energy E from xfLi 
to xf [reversely for T^(E)}. O wing to time reversal invariance we have Tl(E) = Tf (£)[= T^E)}. For Boltzmann 
statistics (nondegenerate regime), we have f(x;e) = (Airm* / f3h 3 ) exp{— j3[e + E c (x) — Ef(x)]}, and obtain from Eq. 
® 
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= ~ev c N c T(i) 



,/3£ F (*™ i) _ p pE F {x?) 



1,...,N; 



(2) 



here, v c = (27tto*/3) _1 / 2 is the emission velocity, the factor N c = 2(27rm*//3/i 2 ) 3 / 2 is the effective density of states 
at the conduction band edge, E F (x) is the quasi- Fermi level at position x, and /? = l/k^T. The factor T(i) is the 
thermally averaged probability for ballistic transmission across the interval {i}, 



/>oo 

T(i)=(3 dE e~P E Ti(E) , 



(3) 



where = max{£? c (x i _ 1 ),^ c (i i )} (cf. Ref. ||, and Ref. ||, Sec. 2.1). 

In the classical description, we have Ti(E) = ®{E — E™(i)), where E™{i) is the maximum of the conduction band 
edge profile -E c (aO between or at the points and asf- Therefore, Eq. (||) becomes 



i(0 



e /3£ F K™ i) _ e /3B F (xf ) 



»=1,...,JV 



(4) 



In what follows, we adhere to the classical description since it allows the greatest transparency, in particular as far as 
the averaging over ballistic intervals is concerned. The full formulation including tunneling effects will be presented 
in Sec. III.B. 

In the stationary case, the current is independent of position, = j = const., and we can write Eq. (|]) in the 
form 



ev Q N c 

Iterating this relation, we find for a configuration with N ballistic intervals 

JV 



1,...,N 



e /3E F (S) _ e {SE F (0) 



J 



,/3S»(i) 



(5) 



(6) 



The sum on the right-hand side must be averaged over all possible configurations of the ballistic intervals, i.e., over 
all positions of their end points (i — 1, TV — 1) in the chain of N intervals, where N = 1, oo. Introducing 
the absolute maximum E™ of the band edge profile E c (x) in the interval [0, S], we denote the product of exp(—(3E™) 
with the average of the sum by H, 



N 



(7) 



\i=l 



Setting E F (0) - E F (S) = eV, we then have 



1 



-f3eV _ 



ev c N c 



(8) 



where E p = E" 1 — E F (0) is the overall barrier height. The end points of the sample are connected to large ideal 
reservoirs in equilibrium characterized by the Fermi levels -Ef(O) and E F (S) whose difference determines the voltage 
bias V \M. 



B. Averaging over the ballistic configurations 

1. The distribution of the ballistic intervals 

The distribution of a set of N ballistic intervals, i.e., of the points of equilibrium xf , across the sample of length 
S is determined in the following way. The probability of an electron to make a collision in the interval g?£ after 
having traversed a distance £ since its last collision is given by exp(— £/l)(d£/l), where I is the mean free path p^| , ^6[ . 
Assuming (as in the Drude model) that an electron which collides with an impurity is taken out of the ballistic 
current and can be counted as equilibrized, we write the complete distribution of the ballistic intervals in the form of 
an infinite-dimensional diagonal matrix with elements labelled by the number N of ballistic intervals, N = 1, 2, ... : 
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dp 



-S/l 



N-l 

n 



dx 



N 



-(S-x%_ 1 )/l 



{...,dP N ,...} ; 



(9) 



here the exponents in the general term cancel out except for the term (—S/l), and we have (recalling Xq = 0) 



W " X dx? 



dP x = e- s ' 1 ; dP N = e- S ' 1 J] ^- 0(xf - xfL,) for JV > 2 



(10) 



This distribution is normalized to unity, since 



Tr [dP=(l+jr i [ 

J \ N=2 Jo 



dr N 



-S/l 



(11) 



Here, the infinite sum over N contains power terms which add up to the exponential exp(S/l); the evaluation of all 
quantities considered in the following runs along similar lines. 

The average sum (Q) can be written as an expectation value of the form 



\i=l / {N,xf} 

where we have introduced the infinite-dimensional diagonal matrix 

( N \ 



(12) 



N = l,2,.. 



(13) 



2. The average for a peak in E c (x) 



The case where the conduction band edge profile E c (x) has a single peak at the position X somewhere along the 
sample is illustrated in Fig. n\ If a ballistic interval {i} contains X, we have E™(i) — E C (X); if it lies to the left [right] 
of X, we have E™(i) = E c (xf l ) [= E c (xfLx)]. Thus we obtain for the average sum 



AT=1 



N 



-0E" 



,/So(xf ) 9{x _ X N ) + e 0E c (X) d{x N _ x) d{x _ ^ + g/JBcCxilO Q{x N^ _ x) 



= l + S/l, 

where E™ = E C (X) in the "reduced sample length" 



(14) 



S= I dxe-K E ?- E °W 
'o 



(15) 



satisfying S < S. Special cases are X = and X = S, when the profile is monotonic. 

It is seen that the average sum S for a profile containing a single maximum at a position X inside or at the end 
of the sample, is given by unity plus the ratio of reduced sample length and mean free path. By the definition of the 
average sum E, the unit term represents the contribution of the ballistic transmission across the highest peak E™ of 
the profile E c {x). 



3. The average for a valley in E c (x) 

We consider a conduction band edge profile of the type shown in Fig. || It contains two peaks at Xq and X\, 
respectively (without loss of generality, the left peak is assumed to be the higher one), and in between a valley with 
minimum at Y\. The average sum over the ballistic intervals enclosed by the peaks (Xq < xf_^\ x? <X\) is given by 
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N=l 



N 

E« 

i=l 



cLPn ? « 



(16) 



where x* is that position to the right (left) of Y\ where the profile has the same height as at the point x to the left 
(right) of Yi. The average sum over the ballistic intervals contributed by the left-hand peak at X is given by Eq. (|l4| ) 
with X = Xq and the last term in the brackets omitted, and analogously for the right-hand peak at X\. The average 
sum for the profile of Fig. || can then be evaluated as 



S/l 



^L e -\x~x'\/l e -l3E s (X ) 



a l3E c (x) 



(17) 



where now E™ = E C (X ) in S. 




FIG. 2. Averaging over the ballistic configurations: single valley of the band edge profile E c (x). For explanation, see text. 



It is seen from this formula that for &o,i C'Ctci, where wi = X\ — X^ is the width of the valley and 6o,i are the 
"widths" of the two barriers at Xo,i, the integral over the second exponential in the brackets can be neglected. Since 

/ 1 dxe-\ x - x "U l = l{l-e- w ^ l )ttl , (18) 
Jxi 

we obtain 

Z = l + S/l + e-^ x °)- E °™ , (19) 

i.e., we have independent contributions from each peak (first and third terms). 

On the other hand, for S <C I (ballistic regime) , the whole integral term in Eq. ( |l7| ) can be neglected (together with 
the term S/l) since it is smaller than (wi/l) exp{— (3[E C (X$) — E c (Xi)]}. Then S reduces to the unit term, which 
represents the contribution of the transmission across the higher peak at Xq: this peak "eclipses" the lower peak at 
Xl 



III. CURRENT- VOLTAGE CHARACTERISTIC 

A. Classical current-voltage characteristic 

In the foregoing, we have considered single peaks and valleys. If the profile contains not just one, but an arbitrary 
combination of such structures as in a chain of grains in polycrystalline materials, the average sum S for a general 
band edge profile is given by 

H = l + (5 + A)//, (20) 
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where we have introduced the "shape term" A for M valleys, 



e 0E c (X v ) _ £ l3E c (x) 



(21) 



The contribution of valley v (with adjoining lower maximum X v ) consists of an integral which extends over the width 
of the valley from X* to X v . We now define the effective transport length L — El, so that 

L = l + S + X. (22) 

In compliance with Eq. (p"8|), we find A < S. 

From Eqs. (0), ( ^o| ) and we obtain the classical current- voltage characteristic for nondegenerate systems as 

j = -ev e N c e-? E * j- (1 - e~^ v ) . (23) 

This formula is the principal result of the present work. It is given here in the context of classical transport, where 
its interpretation is most perspicuous; corrections due to tunneling will be introduced below. 

The properties of the current- voltage characteristic ( p3| ) are determined by the barrier height E p and the ratio l/L. 
In Eq. (|2^) for the effective transport length L, the mean free path I represents the ballistic contribution to the current 
(this term is associated with the highest peak of the band edge profile). The remaining terms give a quantitative 
measure of the influence of that part of the electron motion which is not purely ballistic. Their contribution amounts 
to at most twice the length S of the sample. The reduced sample length S given by Eq. (^|) represents a contribution 
that characterizes the band edge profile E c (x) in an integral way; it does not manifestly depend on I, only indirectly so 
via the profile (cf. below) . The shape term A given by Eq. ( |2~l| ) , on the other hand, depends on the detailed structure 
of the profile as well as explicitly on the mean free path and thus represents the interplay of ballistic and diffusive 
transport. This term is a distinctive feature of formula ([23]). 

We emphasize that the integrals appearing in the effective transport length L, and in particular the reduced sample 
length S, result from averaging over ballistic configurations, not because I were so small that the sum over intervals 
could be replaced with an integral, as assumed in the diffusive regime. 

The barrier height E p and the effective transport length L are determined by the band edge profile E c (x) which is 
a solution of the Poisson equation 

E'>(x) = -[-en(x)+Q(x)}, (24) 

where 

n(x) = Nce -l3[E^)-E F {x)] (25 ) 

is the conduction electron density and Q(x) is the density of fixed charges. In the "trapping model" for grain 
boundaries in polycrystalline materials |24|-|27|] , which our calculations in Sec. IV are based upon, the density Q{x) is 
given by 

M 

Q(x) =eN doa + J2<lU(x-X v ) , (26) 

v=l 

where N^on is the density of donor atoms (assumed completely ionized), and g* is the area density of the charge 
associated with occupied acceptor-like "trapping states" localized at the grain boundary at X v , The donor density 
Ndon not only affects the band edge profile E c (x), but also determines the magnitude of the mean free path I (a 
simple relation between I and Adon is obtained from the analytical expression for the A^don-dependence of the electron 
mobility [i given in Ref. ps{ ], using /i = ev c [3l). Thus, there is, in general, an indirect relation between E c (x) and I, 
which implies an indirect (and generally rather strong) ^-dependence of the barrier height E p as well as of the reduced 
sample length S and the shape term A [as noted above, the latter term also depends explicitly on I, namely via the 
weight function exp(— \x — x*\/l)]. 

When the bias vanishes, the electron density n(x) has the equilibrium form given by expression ( p5[ ) with E-p(x) = 
const. = Ep(0), and we obtain a nonlinear differential equation for E c (x). In the presence of bias, the Poisson equation 
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( p4| ) must be solved in conjunction with an equation for the current j. In the diffusive limit, this equation is given by 
the familiar drift-diffusion expression for the current, 

j = fj,n(x)-^E F (x) = ^^ n ( x ) + fJ-n(x)-^E c (x) , (27) 

which determines n(x) in terms of E c (x) and the constant parameter j. On the other hand, when the current is 
ballistic in an interval (xfLi,xf) [cf. Eqs. (EJ)], only the equilibrium densities n(xf!_i), n(xf) at the end points of the 
interval enter into the description, and an averaging formalism must be provided which allows one to derive from 
these discrete values of the density a continuous physical density to be used in the Poisson equation. 

We close this section with a brief discussion of the special case of a single barrier. In the case of a single grain 
boundary or a Schottky contact, the band edge profile exhibits a single peak and no valleys, so that the shape term 



vanishes, A = 0, and Eqs. (22) and (|23|) yield 



j = -ev c N c e - pE »— L=(l - e~ PeV ) . (28) 
l + S 

For a grain boundary, the barrier height E p is given by the difference of the profile maximum at the boundary and 
the Fermi level in the bulk of the grain. The barrier height of a Schottky contact is equal to the difference of the 
profile maximum in the semiconductor and the Fermi level of the metal. Equation (|2^) is formally identical to Eq. 
(6) of Ref. Q (except for the tunneling correction). However, the present derivation is different from that of Ref. ||, 
the averaging over ballistic configurations being the crucial ingredient. 

If the mean free path is much longer than the reduced sample length, I ^> S, one obtains from Eq. ( p8[ ) the 
thermionic-emission formula EL^. In the opposite case, / -C S, the transport mechanism is diffusive; with the use of 



Eq. (p8|) becomes 



li = ev c pi , (29) 



(3S 

Looking at the diffusive limit differentially, we have from Eq. (Eh in conjunction with Eq. (|25|), replacing the discrete 
coordinate xf with the continuous coordinate x, i.e., setting xf = il = x (i = 1, . . . , N), I — S/N with N — > oo, 

j = pev c N c e^> 3[E ^ x) ~ EF ^ ] [dE F {x)/dx} ■ I = l in{x)^E F {x) , (31) 

ax 



in agreement with Eq. (12 



B. Quantum effects 

Since we are focussing attention on the case where the mean free path and the relevant structures of the band edge 
profile are of comparable length, we are generally dealing with systems of small dimensions and therefore must expect 
quantum effects such as discretization of energy states and tunneling to play a role. Considering, e.g., a potential 
valley with a width of the order of 20 nm, as typified by the examples discussed below, we find that the electron 
motion is quantized with energy spacings of about 0.03 eV. Since the barrier heights in the examples exceed this 
value by an order of magnitude, one may still speak of a classical continuum of states. On the other hand, the wave 
length of the electrons at T = 300 K is A = h/(m*v c ) 10 nm, so that quantum tunneling should be important. The 
formalism for including the effect of tunneling, i.e., for going from thermionic emission to thermionic field emission, 
will be developed in the following. We do not include the effects of phase interference and localization, since these 
are not expected to play an important role for the polycrystalline materials we are considering. 

The transmission probabilities for ballistic intervals with no peaks inside are treated classically, as before. Tunneling 
has to be taken into account near each peak in intervals X~ < x < containing the peak at X n (n = 0, 1, . . . , M). 
Since we can treat only ballistic quantum transport, these intervals must not contain equilibration points. In other 
words, in these intervals we must effectively set I —> oo, and in the integrals over x these intervals are omitted, since 
here dx/l — 0. The lengths to be chosen for the intervals will be discussed below. 
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In WKB approximation, the thermally averaged quantal probability T WKB (i) for ballistic transmission from xf i _ 1 
to a;^ is given by 

T WKB {i) = e ^™« + (3 W dE exp ^-0E - | jT' dx{2m* [E c (x) - £]} 1/2 ^J , (32) 

where E^ is defined after Eq. ([|). The limits of integration yi and y* are the turning points at energy E on either 
side of the peak (if the interval {i} = [xf_ 1 ,xf i ] contains several peaks, the integral in the exponential goes from the 
left-most turning point to the right-most). 

It is found that in the current-voltage characteristic (|j|), the barrier factor exp(— (3E P ) is to be multiplied by the 
tunneling correction C , with 



twkb (x - )X + ) rwKB (x - )Xo+) y rwKB (J -xj 



(33) 



where the notation T WKB (X^ , X^- ) refers to the averaged WKB transmission probability T WKB (i) for the interval 
[Xq,Xjj\i etc; here, Xq is the position of the overall maximum peak E™. The reduced sample length S and the 
shape term A that enter the effective transport length L [cf. Eq. ([2^) 1 now appear as 

S = C- 1 { dxe~^- E ^ , (34) 



i) 



where the upper limit S* implies that all intervals [X n , X£] (n — 0,1, ... , M) are to be omitted, and 

M r x~ 

I = C~ 1 J2 / (fce-l'-'l/'e-^ke^^'-e^'W , (35) 



jj=i 



with X v ' = max{X v and 

A v = T(X-,X+)/T WKB (X-,X+) , v = 1, . . . , M , (36) 

where T is the average classical transmission probability (||). It can be shown that A < S, as in the classical case. 
For a single barrier, Eq. J2S|) is generalized to 

j^-eweJVcC^e-^-Lfl-e-"^) . (37) 

According to Eq. (|34|), 5 includes the tunneling correction factor C _1 > 1, and thus tunneling enhances the relative 
effect of the non-ballistic part of the electron motion embodied in the term S in the denominator of expression (|37|). 

The intervals [X~ , A+] enclosing the peaks at X n are defined as the intervals in which tunneling plays a role; they 
are determined by the requirement that the ratio A v becomes virtually independent of the interval if that extends 
beyond the end points X^. On the other hand, it must be ascertained that these lengths are smaller than the mean 
free path I, otherwise the quantum correction scheme breaks down: if tunneling takes place over distances larger than 
the mean free path, thermalization and tunneling occur simultaneously, which cannot be described in the present 
framework. Explicit calculations of the quantum corrections will be carried out below. 



C. The degenerate case 

In the degenerate case, a simple treatment is possible only in the limit of zero bias. Formula (|2|) for the current is 
to be replaced (for infinitesimal bias) with 

where we have set Ef(x^_ 1 ) = E-p and Ep(xf) = Ep + Si. Classically, this becomes 
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(3h 3 



By analogy with Eq. (o), we then find 



MS) - MO) = £ * = i 4^: E i n( i + e /^ K -^(i)] ) 



(39) 



(40) 



with £^f(0) and E F (S) differing infinitesimally from the equilibrium value E F . Denoting here the average over the 
sum (cf. Sec. II.A), multiplied by ln(l + e ^ EF ~ E ^), by S d , 



I" ln(l + e^^- g °]) 1 



we have 



eV = E F (S) - E F (0) = -j 



Ph 3 



{AT,xf } 



1 



(41) 



(42) 



4nem* ln(l + e^ E *- E ^) 

The calculation of the average sum S d follows Sec. II. B. A (classical) effective transport length L d is introduced as 

L d = S d I = I + S d + A d , (43) 

where 



SV1 = dx 



l D (l + e /3^-g]) 
l n (l + e /3[£ F -Bo(x)]) 

is the reduced sample length for the degenerate case, in generalization of Eq. (|l5|), and 

ln(l + e^- E ^) ln(l + e^- E ^) 



M f x v 
A d =E / d"- 11 - 1 " 1 ^ 



a: ; 



ln(l + ePlEF-E c (x v )]) i n (! + e p[E F -E c {x)\} 



is the generalization of the shape term (|2l|) . The zero-bias conductance per unit area, finally, is obtained as 



= ^fHl + e-^)± 



(44) 



(45) 



(46) 



with E p = E™ 1 — Ep (note that in the degenerate case, E p may be negative). Quantum effects may be taken into 
account as in the nondegenerate case. 



D. The effect of dimensionality 

The samples we consider are assumed to have a parallel-plane structure and thus have an essentially one-dimensional 
geometry. Nevertheless, the treatment of the transport process itself should, of course, be three-dimensional. There- 
fore, in the discussion preceding Eq. (^|), one must consider the distance between collisions, £, in three-dimensional 
space. This introduces a polar angle to be summed over appropriately in averaging over the ballistic configurations, 
thereby rendering the formalism appreciably more complicated. In our effort to obtain a simple, closed formula for 
the current-voltage characteristic, we have restricted ourselves to a one-dimensional formulation, which leads to the 
desired result in a straightforward way. 

Now, for the special case of a constant profile, E c (x) = 0, de Jong j2^| has derived a formula which describes 
the transition from the (ballistic) Sharvin resistance ]3(| to the (diffusive) Drude zero-bias conductance |^5|,|lJ. In 
this work, use is made of an integro-differential equation in two and three dimensions, leading to a numerical result 
which, for the three-dimensional case, is summarized in the formula (in the notation of our Eq. ( ^6|) with E™ = and 
PE F > 1) 
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47re 2 m*i?F 

~h? I + hS 



3= T^T-T-^' (47) 



where 7 increases monotonically from 1 to 4/3 as l/S goes from (diffusive limit) to 00 (ballistic limit). In our 
one-dimensional formulation, the denominator in formula (^?J) reads simply I + S. We therefore suggest that the 
restriction to a one-dimensional formulation introduces a numerical error of perhaps 30 percent; the main features 
of our model, which the following examples show to involve enhancements by factors of 10, should not be severely 
affected. 

Unfortunately, one cannot make use of differential equations in the more general case when the profile E c (x) 
changes significantly over distances of the order of the mean free path I. Anyhow, we are not so much interested in the 
numerical solution of the transport problem (which would most efficiently be handled by Monte Carlo simulations) 
but, as far as possible, in its physical analysis via simple and transparent closed formulas. 



IV. EXAMPLES 



A. Zero-bias conductivity of a chain of identical grains 

The zero-bias conductance per unit area, g = (\j\/V)v^o> is obtained as 

g = (3e\ c N c l - = ^f- l - . (48) 

The band edge profile E c (x) for zero bias is given by its equilibrium shape, and is calculated as a solution of the 
Poisson equation ( p4| ) in the trapping model. 

We consider a chain of v identical grains, each of length s, as shown in Fig. ||. There are v — 1 identical valleys in 
E c (x). Writing S = vs and A = (v — 1)1, we have from Eq. ( ^2|) 

L = I + us + (v - 1)1. (49) 

As mentioned above, the profile E c (x), and therefore also the barrier height E v , the reduced grain length and the 
single-grain shape term I, indirectly depend on the mean free path I via its connection with the donor density A^on- 
The /-dependence of the main feature of the profile, the barrier height E p , can be read from the /-dependence of the 
conductance g for a single grain in the ballistic regime, when L = I in Eq. (f4g|). 




FIG. 3. Equilibrium band edge profile E c (x) of a chain of v identical grains of length s, forming a sample of total length 
5 = vs. 



The term I in expression ( |49| ) for the effective transport length L represents the ballistic contribution to the current. 
For purely ballistic motion (/ vs and L — > I), one finds from Eq. (S) in the classical limit (C = 1) that only one 
barrier, the front barrier, is relevant since it eclipses all others (cf. above, and also the discussion in Sec. II of Ref. 
p2[). However, when tunneling takes place (possibly through several neighboring peaks), Eq. ( |37| ) in conjunction with 
Eq. ( |33| ) is to be used to obtain the conductance. 

For the purpose of illustration, we consider two silicon samples with grain length s = 30 nm (/ic-Si) and s — 100 
nm (pc-Si), respectively. Instead of the zero-bias conductance g, we introduce here the zero-bias conductivity a as 
conductance times sample length, 
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a = gS = gvs 



(50) 



(in the diffusive limit, a is independent of the number of grains v). The conductivity has been calculated at temperature 
T = 300 K by means of formulas (Q) and ( [l9|) as a function of the donor density iVdon- Applying the criteria discussed 
in the final paragraph of Sec. III.B, the length of the "tunneling intervals" [X~,X+] enclosing the peaks was taken 
to be 6 nm (which, not surprisingly, is of the order of the wave length A = 10 nm of the electrons); this value is much 
smaller than the magnitude of the mean free path in all cases considered. The results are shown in Fig. |as a function 
of the mean free path I, employing the connection between I and iVdon obtained from Ref. ]2S|] and displayed in the 
inset. It is seen that the conductivity depends appreciably on the number of grains in the case of the smaller grain 
size, s — 30 nm. Thus it emerges that the conductivity calculated for a single grain cannot simply be carried over to 
the entire microcrystalline sample. The latter must be considered as a whole; there is in this case no grain-specific 
conductivity. 




FIG. 4. Zero-bias conductivity a in the unified model for Si samples with grain length s = 30 nm (upper panel) and grain 
length s = 100 nm (lower panel), plotted as a function of mean free path I for temperature T = 300 K, trapping state density 
Nt = 2 x 10 12 cm -2 , and a single trapping level located at 0.56 eV above the valence band edge. Dashed curves: single grain 
[y = 1); solid curves: chain of a hundred grains (y = 100). The inset in the upper panel shows the relation between I and the 
donor density iVdon at T = 300 K. 



Tunneling has a direct influence on the barrier transmission probability, which is taken into account through the 
correction factor C _1 in Eq. (Esb. It enhances the transmission probability by up to 50% in the region of small I. 

In Fig. |^, we summarize the results of the unified model in comparison to those of the drift-diffusion and thermionic- 
emission models. Since tunneling does not affect appreciably the mutual relation between the different transport 
models, we here consider only the classical conductivities, and plot these relative to the conductivity within the drift- 
diffusion model, odd [with L = us in Eq. (fig)]. The curves labelled UM(1) and UM(100) represent the conductivity 
of the unified model for a chain of one and a hundred grains, respectively, divided by the conductivity calculated 
within the drift-diffusion model. In the regions where these curves approach unity, the transport mechanism is 
predominantly diffusive. The curve labelled TE ("thermionic emission") represents the relative conductivity for 
purely ballistic transport across a single grain boundary, the result being identified (as is usually done) with the 
conductivity of the entire chain. This procedure ignores the eclipsing effect alluded to above, and is justified only if 
the mean free path is long compared to the width of the barrier but short compared to the length of the grain [cf. 
the conditions leading to Eq. (|l9|)], so that while moving through the grain, the electrons are thermalized and "face" 
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each grain boundary barrier with the same thermal distribution as when passing over the previous one. 

We see that the transport mechanism tends to become diffusive for very small values of I, and also for I > 80 nm 
in the case v = 100; this holds for s = 30 nm as well as for s = 100 nm. The thermionic-emission model yields 
acceptable results in the region 20 < I < 60 nm at s = 100 nm for v = 1 and v = 100. At s = 30 nm, it is valid for 
v = 1 but completely invalid for v = 100. In the latter case, neither the drift-diffusion nor the thermionic emission 
models describe adequately the correct mechanism [which is represented by the curve labelled UM(100)]. 
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FIG. 5. Relative conductivities ct/odd for the cases of Fig. 



The transport properties have been discussed here for just two grain lengths s representative of microcrystalline 
and polycrystalline silicon, respectively. A more comprehensive study would have to deal with the s-dependence over 
a suitably broad range, which in particular may lead to the disclosure of possible scaling properties. 



B. Degenerate case: Single barrier 



We apply formula ( 46 ) to transport through a single grain boundary. Choosing parameter values so as to reproduce 
approximately the conditions of the example of Ref . [ pd] | , we consider the temperature dependence of the conductance 
for a grain boundary barrier in highly doped pc-Sn02:Sb with grain length s — 50 nm. We have calculated the 
equilibrium band edge profile E c {x) at T = 300 K and A^on = 7.2 x 10 18 cm~ 3 , corresponding to a mean free path 
I = 11 nm. We assume a single trapping level at midgap (1.75 eV above the valance band edge) with trapping state 
density At = 1.95 x 10 12 cm~ 2 . The barrier height is obtained as E p = 0.038 eV and the barrier width as w 6 nm. 
Since we intend to apply formula ( f46|) in a schematic manner only and the band edge profile is found to depend only 
weakly on T and A*don (or Ep), we adopt the profile calculated at T — 300 K for all values of T and control the degree 
of degeneracy by independently choosing the value of E-p in Eq. (46). 

In Fig. ^, we display the T-dependence of the barrier conductivity a = gS, calculated from Eq. (^) using for the 
mean free path the fixed value I = 11 nm (when comparing with the results of Ref. one must take account of 
the factor a c s « 10~ 2 introduced there). The case considered here [Ep = 3.52 eV, compared to E c (grain bulk) = 
3.50 eV] is strongly degenerate; nevertheless, the results of calculations using Boltzmann and Fermi-Dirac statistics 
[UM(non-d) vs. UM(d)] are not too far apart. Further, it is observed, by comparing the curves labelled UM(d) and 
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UM(d,cl), that the effect of tunneling increases the conductivity by more than a factor of 2, becoming stronger as the 
temperature decreases. 

The purely diffusive and ballistic (degenerate) conductivities are also shown in Fig. || It appears that, as the 
temperature rises, the transport character in the unified model changes from ballistic to diffusive. This is governed by 
the relative size of the terms I and in the effective transport length Ld appearing in the denominator of expression 
( [4(f ) for the conductance g. 

200 
150 



E 

50 


100 150 200 250 300 

T [K] 

FIG. 6. Temperature dependence of the barrier conductivity a for a grain boundary barrier in pc-Sn02:Sb with conduction 
band edge profile calculated for the fixed parameter values s — 50 nm, T = 300 K, JVdon = 7.2 x 10 18 cm' 3 (I — 11 nm), 
Nt = 1.95 x 10 12 cm -2 , and a single trapping level located at 1.75 eV above the valence band edge. The Fermi level is chosen 
as Ef = 3.52 eV. UM(d): unified model for the degenerate case (Fermi-Dirac statistics); UM(non-d): unified model for the 
nondegenerate case (Boltzmann statistics); DD: drift-diffusion model; TE: ballistic model; UM(d, cl): unified model for the 
degenerate classical case (no tunneling). 




Formula (44) yields for a value which increases with Ep, i.e., with donor density N^ on . In Ref. [^T), the 
corresponding term jw is found (by fitting to data) to decrease instead; a more detailed analysis appears to be 
necessary for an explanation of this discrepancy. 



V. SUMMARY AND CONCLUSIONS 



For electron transport in parallel-plane semiconducting structures, we have developed a generalized Drude model 
which unifies ballistic and diffusive transport for arbitrary magnitude of the mean free path and arbitrary shape of 
the conduction band edge profile. The semiclassical approach has been adopted, but tunneling has been taken into 
account in WKB approximation. 

The basic assumption of the model is that the electrons move ballistically over intervals whose lengths are randomly 
distributed about the value of the mean free path. By averaging over the random configurations of ballistic intervals, 
we have derived simple formulas for the current-voltage characteristic (in the nondegenerate case) and for the zero- 
bias conductance (in the degenerate case). The distinctive feature of these formulas is the presence of an effective 
length that comprises a shape term directly manifesting the interplay of ballistic and diffusive transport. Previously 
obtained formulas for the current-voltage characteristic and for the zero-bias conductance refer to special cases and 
do not include such a term. 

We have performed numerical calculations of the zero-bias conductivity for chains of grains of (nondegenerate) 
/iC-Si and pc-Si, and for a single grain boundary in highly doped (degenerate) pc-Sn02:Sb. The calculations for Si 
show substantial deviations of the results of the unified model from those of the purely ballistic and purely diffusive 
models. Moreover, within the unified model, one finds a fairly strong dependence of the conductivity on the number 
of grains. 

For the calculation of the zero-bias conductivity, the band edge profile to be used is the equilibrium profile. Except 
in this case and in the diffusive limit, the determination of the band edge profile is highly complicated in general. 
Basically, one has to solve the Poisson equation along with the relevant current equation self-consistently, taking into 
account the averaging of the electron density over configurations of ballistic intervals. 
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In this work, we have considered transport in one dimension; it has been argued that this should be sufficient to 
exhibit the essential features of the physical phenomena involved. Possible generalizations would be the inclusion of 
the coupling to minority carriers, of optically induced carrier generation, and of recombination. 
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